Towards an Explanation of the Mesoscopic Double-Slit Experiment, 
a new model for charging of a Quantum Dot 
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For a quantum dot (QD) in the intermediate regime be- 
tween integrable and fully chaotic, the widths of single- 
particle levels naturally differ by orders of magnitude. In 
particular, the width of one strongly coupled level may be 
larger than the spacing between other, very narrow, levels. In 
this case many consecutive Coulomb blockade peaks are due 
to occupation of the same broad level. Between the peaks 
the electron jumps from this level to one of the narrow levels 
and the transmission through the dot at the next resonance 
essentially repeats that at the previous one. This offers a 
natural explanation to the recently observed behavior of the 
transmission phase in an interferometer with a QD. 

PACS numbers: 73.23.Hk, 05.45.-a, 73.20.Dx 



In spite of much progress in the fabrication and exper- 
imental investigation of ultrasmall few-electron devices - 
such as Quantum dots |Q , many experimentally observed 
features of these systems still remain unexplained. 

A challenging problem which has resisted adequate 
theoretical interpretation arises from the experiment j^] 
which determines the phase of the wave transmitted 
through the QD §. The main goal of this paper will 
be to find a mechanism which may lead to a satisfactory 
explanation of these results. Hopefully our approach will 
also allow to shed some light on other open problems 
concerning the Coulomb blockade (CB) Q in QD-s. 

In the experiment of ref. in addition to the con- 
ductance of the QD, the phase of the electron transmit- 
ted through the QD was measured via an interference 
arrangement. In accordance with the Breit-Wigner pic- 
ture, the phase increased by tt around each CB peak. 
Absolutely unexpected, however, was a fast jump of the 
phase by —n between the resonances near the minimum 
of the transmitted current. Such a behavior is in evident 
contradiction with what one would expect if the trans- 
mission of the current proceeds via consecutive levels in 
1-dimensional quantum well. 

In a two-dimensional QD the phase drops associated 
with the nodes of the transmission amplitude arise al- 
ready within the single-particle picture P,^. However, 
in order to have a sequence of such events one should 
consider a QD of a very special form. The model of 
ref. 0] also does not allow to explain the series of drops. 
The mechanism of refs. |^,^) makes nontrivial assump- 
tions on the geometry of the QD and the way it changes 



under the change of plunger gate voltage. An interest- 



ing generic mechanism suggested recently in ref. 10 may 
indeed lead to the correlations in transmission at many 
consecutive valleys, but the predicted phase behavior dif- 
fers from what has been seen experimentally. 

In this paper we propose a mechanism according to 
which the transmission at many CB peaks proceeds 
through one and the same level in the QD. This means 
that the phases at the wings of different resonances 
should coincide and the increase by tt at the resonance 
must be compensated. This compensation occurs via 
narrow jumps between the resonances and is accompa- 
nied by a fast rearrangement of the electrons in the dot. 

Although the experiment Q was clearly done in the 
CB regime, the widths of the resonances turned out to 
have been anomalously large, only few times smaller than 
the charging energy. Also the widths and heights of all 
observed resonances are very similar. These surprising 
features of the results of ref. Q , which have not attracted 
so wide an attention as the phase jumps, also find natural 
explanation within our picture. Our mechanism requires 
the QD not to be fully chaotic (neither do we require 
an integrable QD). It is not clear, how chaotic was the 
dot used in the experiment. However, the QD containing 
~ 200 electrons was ~ 50 times smaller than the nominal 
elastic mean free path. Thus, disorder should not have 
been essential for the dynamics of the electrons. 

It is generally believed that the CB is observed only 
if the widths of resonances are small compared to the 
single-particle level spacing in the dot A. This condition 
assumes that couplings of all levels to the leads are of 
the same order of magnitude. However, as we will show 
at the end of this paper, even for nonintcgrablc ballistic 
QD-s the widths of the resonances may vary by orders of 
magnitude. In this case it does not make sense to com- 
pare the width of few broad resonances with the level 
spacing, determined by the majority of narrow, practi- 
cally decoupled, levels. 

A useful theoretical model for the description of charg- 
ing effects in QD-s is the tunneling Hamiltonian in the 
constant interaction (Ucb) approximation (see e.g. [H) 
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+ ^2 e(k)a+a k + ^[i;c+a fe + h.c] . 



Here c(c + ) and a(a + ) are the annihilation(crcation) oper- 
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ators for electrons in the dot and in the lead and £j, e(k) 
are the single-particle energies. We do not introduce the 
k dependence of the tunneling matrix elements U. Since 
our approach is mainly based on the energetics of the 
QD, it is enough to consider only one lead. Summation 
over spin orientations is easily included. Also under the 
assumption of capacitive coupling to the gate, the levels 
in the dot flow uniformly with the voltage 

£i = Ei{Vg =0)-Vg. (2) 

The energies of the electrons in the wire are given by 

e(k) = k 2 /2m - E F . (3) 

Here k = nn/L, a (very)large n is the level number in 
the wire and L is the length of the wire. 

For our purposes, it will be possible to simplify further 
the Hamiltonian (Q). We will consider the case where 
the coupling of one particular level N is dominant 
ti, i 7^ N. If the width of this level is larger than the 
single-particle level spacing A, a very nontrivial regime of 
charging of the QD may be described by means of second- 
order perturbation theory estimates. Surprisingly, this 
simple limit of CB has not been considered yet. 

An example of a system for which the widths differ 
drastically is the integrable QD However, it is hard 

to believe that the large (N e ~ 100 1000) QD may be 
even close to integrable. Nevertheless, at least in clas- 
sical mechanics, a considerable gap is left between inte- 
grable and fully chaotic systems. Even in a nonintegrable 
dot two kinds of trajectories - quasi-periodic and chaotic 
may coexist. In this case, in 2-dimensions any trajectory 
(even a chaotic one!) does not cover all the phase space 
allowed by energy conservation. Consequently, the cor- 
responding wave functions do not cover all the area of 
the QD. If such a regime is realized in QD-s, it easily 
explains why the widths of the resonances may vary by 
orders of magnitude. Moreover, many other features of 
such a QD may differ strongly from those of the chaotic 
QD |l7j| . An explicit numerical example, which supports 
the existence of such regime will be given later. 

Now we turn to the many-particle effects arising for 
the Hamiltonian ([!]) in the case of only one (iV-th) level 
in the dot coupled strongly to the wire 



r = T N = 2ir\t N \ 2 dn/de > A 



(4) 



(Here n is the same as in the eq. (g) and dn/de is taken at 
the Fermi energy e = 0.) The widths of the other levels 
are taken to be much smaller than the level spacing and 
may be neglected in the first approximation. The charg- 
ing energy is still very large Ucb 3> T. We shall show 
that transmission of a current at about (r/A) h\(UcB /r) 
consecutive CB peaks will proceed through one and the 
same level e^r. 

Let the levels with i < in the QD be occupied. Our 
aim is to find the total energy E to t of the true ground 



state of the dot at different values of V g . Without loss of 
generality we may assume that the summation over i in 
eq. ([!]) goes only over i > 0. (Thus we subtract from the 
total energy the trivial constant corresponding to selfin- 
teraction of electrons with i < 0. Coulomb interaction 
between electrons at the levels with i < and i > is 
included into £j>o-) Also let us subtract from the to- 
tal energy the trivial energy of electron gas in the leads 
£ £ (fc). 

Let us consider spinless electrons. For large positive 
eN(V g ) ^> T the only contribution to the total energy 
Etot is given by the second order correction (the levels 
in the wire are lowered due to the repulsion from the 
unoccupied level £jv) 
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Here and everywhere below, en (as well as e.g. E\) is 
the function of V g (0). The upperscript (0) at E^ t shows 
the number of electrons at the narrow levels (with i > 0) 
in the QD. Generalization of this result for the case of 
negative En, £jv "C — T (the broad level being below the 
Fermi energy) is straightforward (note that the level En 
is occupied, not the level E\ as one might expect): 



iT,(o) T / 4E F 

Etol = eat - — In j— 
2n \\eN\ 



(0) 



Here the first term £^v accounts for the energy loss due 
to replacement of one electron from lead to the dot. The 
second order level shift now includes both lowering of 
levels with e(k) < en and raising of those with e(k) > £at. 
The perturbative treatment fails for \e(k) — e jv ^ T, but 
the corresponding shifts of levels below and above e^ 
evidently compensate each other, which is equivalent to 
taking the principal value of the integral in Eq ( p|). An 
approach related to ours was used recently in ref. ]ll| for 
the calculation of CB peaks positions. 

Finally, the exact solution (for spinless electrons) for a 
single state interacting with a continuum is also known 
(e.g. |0|). A precise treatment of this situation, along 
the lines of ref. Il5) , yields: 
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which coincides with the Eqs. flqM) at |ejy| 3> T. 

Let us now consider the branch where the level 1 in 
the QD is occupied. The energy of this electron is e\. 
However, adding one more electron via the hopping ijv 
now costs eN + UcB- The ensuing reduction of the down- 
ward shift of the level E^} is of crucial importance. The 
analog of Eq. (pi) for £jy + Ucb > now reads 



E, 



(i) 



•«- Sl 2r: [n \e N + UcB 



(8) 



2 



For small V g one has E^ t < e£} and the Eqs. (|j7fl|) 
describe the true ground state of the system. However, 
the two functions E^l(V g ) and E^ t (y g ) cross at 

UcB 

£N ~ exp{27r( ffJV -£i)/r} + l W 

and the ground state jumps onto the branch E^ t . The 
energy of the current-transmitting virtual state N is pos- 
itive again. Thus, the transmission amplitude phase had 
returned to what it was before the process of filling of 
state N and the subsequent sharp jump into the state 
where level 1 is filled, ft is the latter jump which provides 
the sharp drop by tt of the transmission phase, following 
its increase by tt through the broad resonance. Many 
(~ (T / A)ln(UcB /F)) consecutive resonances are due to 
the transition via one and the same level N. 

For electrons with spin, the Breit-Wigner-related for- 
mula (^) does not work. However, far from the resonance 
the perturbation theory may still be used (at least un- 
til the temperature is high enough to be away from the 
Kondo effect Jig]). We are not able to discuss in detail the 
role of spin in this short note. Still in this case the many 
charging events proceed via one and the same broad res- 
onance, each accompanied by the increase of phase by 7r 
which is compensated by the — ir jump in the valley. 

To illustrate the relevance of the model eq. (Q) with sin- 
gle strongly coupled level we performed numerical simula- 
tions for a model QD of a size I with a simple polynomial 
potential (a smooth QD coupled to two leads) 

V- -4^(1- £)' + (, + £)'(l + 8(f -I)")- (10) 

Due to the strong mixing of the x and y coordinates the 
dot is expected to be nonintegrable, but, similarly to the 
experimental geometry [El, it is approximately symmet- 
ric. For simulations we considered the QD on the lattice 
and used I — 10 which was equivalent to 50 lattice spac- 
ings. The kinetic term is given by the standard nearest 
neighbor hopping. Below we present the results of cal- 
culations with the hopping matrix element r = 18 which 
corresponds to the dot with ~ 100 electrons or ~ 200 
if the spin is included (similar numbers to those in the 
experiment). We have used the potential V of Eq.(|l(i|) 
for < x < I. The lead formed by the potential V = 3y 2 
was attached at x < and a hard wall at x = I. 

Within the energy interval 1.5 < e < 4.7 only one mode 
may propagate along the lead. The analysis of solutions 
of the Schrodinger equation within this interval allowed 
us to find the positions and widths of quasi-stationary 
levels in the dot. As we expected, the widths fluctuate 
very strongly from level to level (by many orders of mag- 
nitude). In particular the widths of two levels #102 and 
#108 exceed sufficiently the level spacing r/A w 6 (the 
number of states doubled due to spin). The widths of 
other levels vary from r/A - 1 to r/A ~ 10~ 5 - 10~ 6 . 



The origin of the hierarchy of widths becomes clear 
from fig. 1, where we have plotted \ip\ 2 in the QD for 
(real) e at the top of corresponding resonances. The 
quantized version of different variants of classical motion 
may be found on this figure. The most narrow level #103 
corresponds to a short stable transverse periodic orbit. 
Other broader levels, such as #96, 106, may be consid- 
ered as the projections of the invariant tori corresponding 
to quasi-periodic classical motion. This classical trajec- 
tory reaches the line V(x, y) — e only at few points. The 
candidates for chaotic classical motion (e.g. #110) also 
correspond to relatively broad resonances |L2]]. Even in 
this case only a part of the QD is covered by the tra- 
jectory. For the most coupled levels #102 and #108 the 
area covered by the trajectory touches the lead by its 
corner. 
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FIG. 1. The density of electrons in the dot at the reso- 
nances coupled to the single-channel wire(attached from the 
left). The line V — is shown dashed. Numbers correspond 
to the number of the level in the QD. About 95% of the norm 
of the wave function in the dot is shown. The "twin copies" 
(such as 100 and 112) of levels 94, 96, 101, 103 are not shown. 

Moreover, two well coupled trajectories contribute to 
the level #102. This is seen from the fig. 2 where we 
show also the \ip\ 2 at the left and right wings of this reso- 
nance. One contribution corresponds to the strongly cou- 
pled quasi-periodic trajectory (left), having the "turning 
point" V(x, y) — e just at the left contact. The other con- 
tribution comes from the true periodic trajectory (right). 
Two quantum states in the dot become mixed via inter- 
action with the wire and form one broad (#102) and one 
almost decoupled (#104) resonance fl3f| . 

We have repeated the calculations several times for 
slightly different V and in a broad range of variation of 
the hopping. Typically we saw the resonances of very 
different width and the origin of the most broad peaks 
was explained by simple classical arguments. 
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FIG. 2. Decomposition of the level 102 into parts corre- 
sponding to simple classical trajectories. Numbers are the 
energies for which the figures were done. 

Taking into account the different sensitivities of lon- 
gitudinal and transverse modes to the plunger |^,|| may 
allow to keep our broad level Spj even longer within the 
relevant strip of energy. This may provide an explanation 
of even longer sequences of resonances accompanied by 
the — 7r jumps. In a more refined approach, adding new 
electrons into the QD should cause a change of the self- 
consistent potential V(x, y). The total energy of the dot 
and the wire will be lowered in the presence of strongly 
coupled levels. This may cause the potential of the QD 
to automatically adjust to allow such levels, which will 
support our explanation of the experiment of Ref . Q . 

Our mechanism of charging of the QD requires the ex- 
istence of the broad level with V ^> A. The simple way 
to justify the relevance of our theory for the explanation 
of the experiment of ref. ^| will be to close the dot suf- 
ficiently in order to have V <C A for all levels. In this 
case the phase still increases by n at any resonance, but 
the correlation between peaks will disappear. (More pre- 
cisely the pairs of peaks corresponding to adding of elec- 
trons with opposite spins onto one and the same level 
still are correlated, but correlation between pairs should 
disappear.) Moreover, within our mechanism a series of 
~ (r/A) \ii(Ucb /F) strong charging peaks in the con- 
ductance should have the same height. This "coupling 
dependent" correlation of the peak heights seems also 
easy to measure. 

To conclude, we have considered the model, for which 
upon increasing V g , it is energetically favorable to first 
populate in the dot the level strongly coupled to the 
leads. At a somewhat larger V g a sharp jump occurs to 
a state where the "next in line" narrow level 1 becomes 
populated. This jump accounts for the sharp decrease 
by ~ 7T of the transmission phase. The similar strengths 
of resonances seen in the experiment g] and their large 
width are also clear within our mechanism. The current 
transmission through such a QD resembles the behavior 
of rare earth elements, whose chemical properties are de- 
termined not by the electrons with highest energy, but 
by the "strongly coupled" valence electrons. The over- 
lapping of single-particle resonances may take place also 
in the Kondo experiments in QD-s Q| , where in order to 
increase the Kondo temperature the dot is usually suffi- 



ciently opened. Hopefully the unusual effects observed in 
some of these experiments may be also explained within 
our approach. 
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